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We use neutron scattering to study the effect of uniaxial pressure on the tetragonal-to- 

orthorhombic structural (Ts) and paramagnetic-to-antiferromagnetic (Tn) phase transitions in 

NaFeAs and compare the outcome with similar measurements on as- grown and annealed BaFe2As2. 

Cn I In previous work on as- grown BaFe2As2, uniaxial pressure necessary to detwin the sample was found 

to induce a significant increase in zero pressure Tn and Ts. However, we find that similar uniaxial 

pressure used to detwin NaFeAs and annealed BaFe2As2 has a very small effect on their Tn and Tg. 

^ I Since transport measurements on these samples still reveal resistivity anisotropy above Tn and Ts, 

f\ , we conclude that such anisotropy cannot be due to uniaxial strain induced Tn and Ts shifts, but 

I ■ must arise from intrinsic electronic anisotropy in these materials. 
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PACS numbers: 74.70.Xa, 74.62.Fj, 75.40.Cx, 75.50.Ee 



^ ■ The parent compounds of iron pnictide superconductors such as NaFeAs and BaFe2As2 exhibit a tetragonal-to- 

orthorhombic lattice distortion at temperature Tg and paramagnetic-to-antiferromagnetic phase transition at Tn 
(^ Ts), forming a low-temperature collinear antiferromagnetic (AF) state with ordering wave vector along the [±1,0] 
' O . directions of the orthorhombic lattice [Figs. 1(a) and 1(b) ] [iHzl- Because of the twinning effect in the orthorhombic 
C ■ AF state, AF Bragg peaks from the twinned domains in Fig. 1(c) should occur at [±1,0] and [0, ±1] positions in 
reciprocal space [Fig. 1(d)] [8|. To probe the possible electronic anisotropic state (the electronic nematic phase) that 
breaks the C4 rotational symmetry of the paramagnetic tetragonal phase in iron pnictides [9| , one needs to prepare 
single domain samples by applying a uniaxial pressure (strain) along one-axis of the orthorhombic lattic e |lOl lll| . 
h^ ■ Indeed, transport measurements on uniaxial pressure detwinned samples of NaFeAs [12| and BaFe2As2 [l3| reveal 
tH- . clear resistivity anisotropy above the zero pressure Tn and Tg that has been interpreted as arising from the spin 
Q> ■ nematic phase |l4l4l6j or orbital ordering |17H23| in the paramagnetic tetragonal state. However, recent neutron 



r^ . scattering experiments on as-grown BaFe2As2 find that a uniaxial pressure necessary to detwin the sample can also 



in 

"^" , the stress-free Tn and Tg in detwinned samples [l2|, ll3| may actually occur in the AF ordered orthorhombic state 



induce a significant (~10 K) upward shift in Tn and Tg l24i] , suggesting that the observed resistivity anisotropy above 



^^ ■ below the strain- induced Tn and Tg. Furthermore, the resistivity anisotropy above Tn and Tg in as-grown BaFe2As2 



and electron-doped BaFe2-a?CoajAs2 becomes much smaller in annealed samples [25|, |26|, suggesting that the observed 
resistivity anisotropy in the tetragonal phase is not intrinsic to these materials but arises from the anisotropic impurity 
scattering of Co-atoms in the FeAs layer [26, 27]. 

Here we use neutron scattering to study the uniaxial pressure effect on magnetic and structural phase transitions 
;J] , in NaFeAs ^], as-grown, and annealed BaFe2As2 ^]. While our measurements on as-grown BaFe2As2 confirm the 
earlier work that the uniaxial pressure necessary to detwin the crystal also causes significant increases in Tn and Tg 
[24] , we find that similar uniaxial pressure has a very small effect on the magnetic and structural phase transitions in 
NaFeAs and annealed BaFe2As2. Since transport measurements on identical NaFeAs and annealed BaFe2As2 show 
clear resistivity anisotropy at temperatures well above the Tn and Tg under uniaxial pressure, we conclude that the 
resistivity anisotropy seen in detwinned NaFeAs and annealed BaFe2As2 in the paramagnetic tetragonal phase must 
be intrinsic properties of these materials, and cannot be due to the impurity scattering effect [26|, |27j]. These results 
suggest the presence of an electronic nematic state in the paramagnetic tetragonal phase unrelated to the Co-impurity 
scattering in electron-doped BaFe2-a?Coa^As2 family of materials [26|, |27|. 

Figures 1(a) and 1(b) show the schematic lattice and magnetic structures of NaFeAs and BaFe2As2, respectively 
[8|. On cooling from the high-temperature tetragonal state, NaFeAs exhibits a tetragonal-to-orthorhombic structural 
transition at T5 ^ 58 K and then a paramagnetic-to-antiferromagnetic transition at Tn ~ 47 K [4]. For comparison, Tg 



and Tat in BaFe2As2 occur almost simultaneously below about 138 K [6|,|7|. In the absence of uniaxial pressure, the low- 
temperature magnetic and crystal structures have equally populated twinned domains with mixed AF orthorhombic 
states as shown in Fig. 1(c). Figure 1(d) shows the [H,K] plane of the reciprocal space where the AF and crystalline 
lattice Bragg peaks for a twinned sample are seen at (±1, 0)/(0, ±1) and (±2, 0)/(0, ±2) positions, respectively. Upon 
applying uniaxial pressure along the orthorhombic ao/bo direction ^], a single domain with sufficient large size can 
be achieved [Fig. 1(e)], the resulting AF Bragg peaks now occurring predominantly at (±1,0) positions [Fig. 1(f)]. 

We prepared high quality single crystals of NaFeAs, as well as as-grown and annealed BaFe2 As2 crystals, using the 
self- flux method. The samples were cut to squared shapes along the ao/bo directions of the orthorhombic structure 
and fit into the aluminum-based detwinning devices for both transport and neutron scattering experiments [Fig. 1(g)]. 
Figure 1(h) shows the comparison of transport measurements for both twinned and detwinned NaFeAs. Consistent 
with earlier measurements ^12,], we find clear resistivity anisotropy below about T* ^ 70 K. Our neutron scattering 
experiments were carried out using the BT-7 triple-axis spectrometer at NIST center for neutron research (NCNR) 
[28], and HB-IA at High-Flux-Isotope-Reactor (HFIR), Oak Ridge National Laboratory. For HB-IA measurements 
on BaFe2As2, the collimations are 48'-48'-sample-40'-68'. The magnetic measurements for NaFeAs were carried 
out on BT-7 with open-50'-sample-50'-120' with Ef = 14.7 meV. To separate the (2, 0, 0)/(0, 2, 0) nuclear Bragg 
peaks in the orthorhombic state of a twinned sample [Fig. 1(d)], we used tight collimation of lO'-lO'-sample-lO'- 
25' on BT-7 with Ef = 14.7 meV. For NaFeAs, the lattice parameters are ao = 5.589, bo = 5.569 and c = 6.991 
A y. BaFe2As2 has lattice parameters ao ^ bo ^ 5.595 A and c = 12.92 A [6]. The wave vector Q in three- 
dimensional reciprocal space in A~^ is defined as Q = i^a* + Kh^ + Lc*, where H^ K^ and L are Miller indices 
and a* = ao27r/ao,b* = bo27r/6o,c* = c27r/c are reciprocal lattice units (rlu). We aligned the crystals in the 
[H, 0, 0] X [0, 0, L] scattering plane, where AF Bragg peaks occur at (±1, 0, L) with L = ±0.5, ±1.5, • • • for NaFeAs Q 
and L = ±1, ±3, • • • for BaFe2As2 [6]. If the twinned domains are equally populated in the zero pressure state, AF 
Bragg peak intensity at (±1,0, L) should be the same as that at (0,±1,L). On the other hand, if uniaxial pressure 
completely detwinnes the sample, the magnetic scattering intensity at (±1, 0, L) in the detwinned state should increase 
by a factor of two compared with the twinned state. 

The detwinning device we used is shown in Fig. 1(g). By knowing the compressibility of the spring and the area of 
the sample, we can estimate the applied uniaxial pressure. For NaFeAs, we have Pq = 0, Pi ~ 7 MPa, and P2 ~ 15 
MPa. The applied uniaxial pressures are Pi ^ 7 and Pi ^ 6 MPa for the as-grown and annealed BaFe2As2 crystals, 
respectively. We first describe the effect of uniaxial pressure on NaFeAs. For our measurements, we always apply the 
pressure at room temperature, and measure the (1,0,1.5) magnetic and (2,0,0) nuclear Bragg peaks. Figure 2(a) 
shows the temperature dependence of the magnetic (1,0, 1.5) peak intensity normalized to the (2,0,0) nuclear peak. 
At zero pressure, we see a clear magnetic intensity increase below Tn = 47 K. On increasing to Pi and then to P2, we 
see that the magnetic scattering intensity almost doubles, suggesting that the uniaxial pressure has indeed detwinned 
the sample. However, the Neel temperatures of the system remain unchanged at Tn = 47 K within the errors of our 
measurements. The normalized magnetic order parameter in Fig. 2(b) shows almost identical behavior for Pq, Pi, and 
P2, thus confirming that the uniaxial pressure needed to detwin NaFeAs has no measurable impact on T^. Figures 
2(c) and 2(d) show wave vector scans along the [i^, 0, 1.5] and [1,0, P] directions at Pq, Pi, and P2. Consistent with 
the order parameter data in Fig. 2(a), the effect of uniaxial pressure is to increase the intensity of the AF Bragg peak 
(1,0, 1.5). To probe the effect of uniaxial pressure on the tetr agonal- to-orthorhombic lattice distortion temperature 
T5, we studied the temperature dependence of the lattice orthorhombicity on the (2, 0, 0)/(0, 2, 0) nuclear Bragg peaks 
using tight collimations. The full-width-half-maximum (FWHM) of (2, 0, 0)/(0, 2, 0) shows a clear increase below 
Tg = 58 K for Pq, Pi, and P2, and thus suggests that the applied uniaxial pressure also has only a small impact on 
Tg [Fig. 2(e)]. Figure 2(f) shows the temperature dependence of the lattice orthorhombicity 6 = {ao — bo)/{ao ± bo) 
at zero pressure and its comparison with the centers of 0-20 scans. We can see a small (~4 K) increase in Tg when 
NaFeAs is detwinned [Fig. 2(f)]. 

Having established that the uniaxial pressure needed to detwin NaFeAs has only a small impact on Tn and T5, we 
investigate the effect of uniaxial pressure on Tn and Tg in as-grown and annealed BaFe2As2. From previous work on 
as-grown BaFe2As2, we know that the uniaxial pressure necessary to detwin the sample will also increase the onset of 
Tn and Tg by ~12 K [24]. On the other hand, transport measurements on as-grown and annealed BaFe2As2 suggest 
that the large resistivity anisotropy in detwinned as-grown samples is due to disorder in these materials and annealing 
significantly reduces the resistivity anisotropy [26] . To determine how uniaxial pressure affects as-grown and annealed 
BaFe2As2, we prepared annealed samples following Ref. [26|, and carried out Tn and Tg measurements on HB-IA. 
Figure 3(a,c) compares the low-temperature normalized AF (1,0,3) Bragg peak intensities at Pq = and Pi = 7 
MPa for the as-grown BaFe2As2. While the overall magnetic intensity behaves similarly with and without uniaxial 
pressure, we see a clear increase in Tn from '^139 K at Pq = to ~141 K at Pi = 7 MPa, consistent with the earlier 
work [2^]. For the annealed BaFe2As2, similar measurements showed almost identical magnetic order parameters [Fig. 



3(b,d)] and a smaller shift in T^ from ^140 K at Pq = to ^141 K at Pi = 6 MPa. Figure 3(e) plots the temperature 
dependence of the FWHM of the nuclear (2, 0, 0)/(0, 2, 0) Bragg peak for the as-grown BaFe2As2. At Pq = 0, the 
peak width increases abruptly below Tg ~ 140 K, reflecting the fact that a twinned orthorhombic crystal has slightly 
different lattice parameters for (2,0,0) and (0,2,0). Pi = 7 MPa uniaxial pressure clearly increases the onset T/v as 
shown in Fig. 3(c), while the Tg under pressure only increases marginally to Tg ~ 141 K [Fig. 3(e)]. Similarly, we 
find that the uniaxial pressure of Pi =6 MPa on annealed BaFe2As2 only increases Tg from ~140 K to ~143 K [Fig. 
3(f)]. 

From the experimental data discussed above, it is clear that the uniaxial pressure necessary to detwin NaFeAs and 
annealed BaFe2As2 has limited impact on T^ and Tg. Theoretically, it has been argued that a small uniaxial strain 
of magnitude Aq should induce an increase in the magnetic ordering temperature ATat = \Aq\ '^ if Tg = Tat, and 
ATat = {Tg-T^Y^ |Ao| if T5 > Tat, where the susceptibility exponent 7 = 2 + 0(l/A/') (with A/" = 3 corresponding to 
the physically relevant Heisenberg case) [29|. The structural transition temperature is also expected to increase on a 
scale of AT5 ~ l^oT, where for A/" ^ 00, x = 1 + 0(1/^") [29|. Comparing with the nearly simultaneous structural and 
magnetic phase transitions in BaFe2As2 [Sl, [31, the structural and magnetic phase transitions in NaFeAs are separated 
by Tg —T/v ~ 12 K [Fig. 1(h)] [4|. Within the spin nematic phase scenario 1 141 . this arises because the c-axis magnetic 
exchange coupling in NaFeAs is much smaller than that of BaFe2As2 [30|, |31|. As a consequence, the shift of the Neel 
temperature ATat = (T^ - Tat)"^ \A^\ ~ \A^\ /144 in NaFeAs should be much smaller than that (ATat = \A^f'^) 
in BaFe2As2, while the changes in structural transition temperatures (AT^'s) should be similar for both materials. 
Indeed, while uniaxial strain seems to have some small effect on Tg for both NaFeAs [Fig. 1(f)] and BaFe2As2 [Figs. 
3(e) and 3(f)], it has virtually no effect on ATat for NaFeAs. This means that the resistivity anisotropy seen in 
NaFeAs above Tat and Tg in Fig. 1(h) cannot be due to the orthorhombic lattice structure or collinear AF order. 
Since NaFeAs also does not have Co as a source for anisotropic impurity scattering [26|, l27| and the Na deficiency 



out of the FeAs plane is not expected to affect the transport measurements [32] , the resistivity anisotropy above Tg 
in NaFeAs must be an intrinsic property of the paramagnetic tetragonal phase under uniaxial strain. These results, 
together with the uniaxial pressure effect on as-grown and annealed BaFe2As2, suggest that the resistivity anisotropy 
in iron pnictide parent compounds cannot arise from strain- induced shift in Tat and Tg. 

In summary, we have shown that the uniaxial pressure needed to detwin NaFeAs and annealed BaFe2As2 has a 
very small effect on magnetic and structural phase transitions of these materials, while transport measurements on 
identical materials reveal clear resistivity anisotropy above Tg. We conclude then that the resistivity anisotropy is an 
intrinsic property in the uniaxial-strained paramagnetic tetragonal phase of NaFeAs and BaFe2As2. 

We thank Jiangping Hu and Weicheng Lv for helpful discussions. The single crystal growth and neutron scattering 
work at UTK is supported by the US DOE BES No. DE-FG02-05ER46202. Work at the lOP, CAS is supported by 
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FIG. 1: (Color online) The lattice and magnetic structures of (a) NaFeAs and (b) BaFe2As2. While BaFe2As2 has the 
orthorhombic lattice and magnetic unit cells, NaFeAs consists of two orthorhombic chemical unit cells stacked along the c-axis. 
The real space schematics of a (c)twinned crystal and a (e) detwinned crystal, the two sets of domains have the same population 
for the twinned crystal whereas one set of domains dominates in the detwinned crystal. In reciprocal space, the magnetic and 
structural peaks corresponding to the two sets of domains have equal intensities for a (d) twinned crystal, while for a (f) 
detwinned crystal the dominant set of domains is enhanced while peaks corresponding to the minority set of domains have 
diminished intensities. Green spheres represent the magnetic (1,0, L) peak and its equivalent points, red spheres represent 
the structural (2,0,0) peak and its equivalent points. The blue arrows in (e) and (f ) represent applied uniaxial pressure, (g) 
Schematic of the pressure device used in this work. Springs of known force constants and area of the sample edge were used to 
estimate the applied pressures, (h) Resistivity of twinned and detwinned NaFeAs. The dashed lines represent Tn determined 
from neutron scattering and T* the onset temperature of resistivity anisotropy, the shaded region represents the temperature 
range of Ts under different pressures found from neutron scattering results. 
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FIG. 2: (Color online) Elastic neutron scattering measurements on NaFeAs under ambient condition and with applied uniaxial 
pressures. Po, Pi and P2 represent ambient condition, ^7 MPa applied uniaxial pressure and ^15 MPa applied uniaxial pressure, 
respectively. The applied pressures are estimated from the changes in length of the spring and the known spring constant, (a) 
Background subtracted magnetic order parameters measured at the (1,0, 1.5) peak normalized to the (2,0,0) structural peak. 
(b) Magnetic order parameters normalized at base temperature show within statistics (uncertainties represent one standard 
deviation) of the current measurement that uniaxial pressure does not affect its shape, (c) Background-subtracted [i/, 0, 1.5] 
scans for (1,0, 1.5) measured at 2.5 K normalized to (2,0,0). (d) Background-subtracted [1,0, L] scans for (1,0, 1.5) measured 
at 2.5 K normalized to (2, 0, 0). (e) Full width half maximum (FWHM) of 0-20 scans at (2, 0, 0) fit with a single Gaussian as a 
function of temperature. The solid lines are guides to the eye and the shaded region is the range of Ts determined from panel 
(f). (f) Orthorhombicity of NaFeAs under ambient conditions determined from fits by two Gaussians of equal intensities and 
centers of 0-20 scans found from fitting a single Gaussian as a function of temperature. The solid lines are guides to the eye 
and the shaded region is the temperature range for Ts. 
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FIG. 3: (Color online) Elastic neutron scattering measurements on as- grown and annealed BaFe2As2 under ambient conditions 
(Po) and with applied uniaxial pressure (Pi is ^7 MPa for the as-grown sample and ~6 MPa for the annealed sample). The 
background subtracted and normalized magnetic order parameters measured at (1,0,3) for (a) as-grown and (b) annealed 
BaFe2As2 under ambient condition and with applied uniaxial pressure, (c) and (d) show expanded plots of the magnetic order 
parameter near the magnetic transition temperature. The arrows indicate temperatures at which the intensity reaches 1% of 
the intensity at 2 K. Full width half maximum(FWHM) of 0-20 scans at (2, 0, 0)/(0, 2, 0) for the (e) as-grown sample and the 
(f) annealed sample fit with single Gaussians as a function of temperature. 



